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ABSTRACT 

The bright Kl III-IV star 7 Cep has been reported previously to have a companion in a ~2.5-yr 
orbit that is possibly substellar, and also has a stellar companion at a larger separation that has never 
been seen. Here we determine for the first time the three-dimensional orbit of the stellar companion 
accounting also for the perturbation from the closer object. We combine new and existing radial 
velocity measurements (of both classical precision and high precision) with intermediate astrometric 
data from the Hipparcos mission (abscissa residuals) as well as ground-based positional observations 
going back more than a century. The orbit of the secondary star is eccentric (e — 0.4085 ± 0.0065) 
and has a period P — 66.8 ± 1.4 yr and a semimajor axis of 19.02 ± 0.64 AU. We establish the 
primary star to be on the first ascent of the giant branch, and to have a mass of 1.18 ± 0.11 Mq, 
an effective temperature of 4800 ± 100 K, and an age around 6.6 Gyr (for an assumed metallicity 
[Fe/H] = -fO.Ol ± 0.05). The unseen secondary star is found to be an M4 dwarf with a mass of 
0.362 ±0.022 Mq, and is expected to be ~8.4 mag fainter than the primary in V and ^--^6.4 mag fainter 
in K. The minimum mass of the putative planetary companion is MpSinz = 1.43 ± 0.13 Mjup, the 
inclination angle of its orbit being unknown. Taking advantage again of the high-precision Hipparcos 
observations we are able to place a dynamical upper limit on this mass of 13.3 Mjup at the 95% 
confidence level, and 16.9 Mjup at the 99.73% (3(t) confidence level, thus confirming that it is indeed 
substellar in nature. The orbit of this object (semimajor axis 1.94 ± 0.06 AU) is only 9.8 times 
smaller than the orbit of the secondary star (the smallest ratio among exoplanet host stars in multiple 
systems), but it is stable if coplanar with the binary. 

Subject headings: binaries: spectroscopic — binaries: visual — planetary systems — stars: individual 
(7 Cep) — stars: late-type 



1. INTRODUCTION 

The bright evolved star 7 Cephei {V = 3.21, SpT 
= Kl III-IV, a = 23'^39"21!01, S = +77°37'55'.'2, 
J2000; also known as HD 222404, HR 8974, HIP 116727) 
is among the first objects to be subjected to high- 
precision radial- velocity measurements in an effort to dis- 
cover s ubstellar-mass companions around nearby stars 
((Campbell. Walker & Yang 1988). This group of inves- 
tigators ( Campbell fc Walker.,.197^ pioneered the use of 
a hydrogen fluoride gas absorption cell on the Canada- 
France-Hawaii telescope (CFHT) and achieved internal 
errors around 13 ms~^ for bright stars, inaugurating the 
era of Doppler searches that has been so successful in 
finding extrasolar planets in the last 10 years. 

Smal l radial-velocity variations in 7 C ep were indeed 
seen bv lCamnbell. Walker fc Yimel l)1988D suggesting the 
presence of a Jupiter-mass object in a '^2.5-yr orbit. 
Those variations with a semi-amplitude of only about 
25 ms~^ were superimposed on a much larger varia- 
tion caused by a previously unnoticed stellar compan- 
ion with a period of decades. However, the interpreta- 
tion of the residual 2.5-yr variation as due to a plane- 
tary object was subsequently put in doubt by th e same 
group (see llrwin et all 119891: iWalker et aP 119891 Il992|) . 
They argued that changes with a similar period were ob- 
served in a chromospheric activity indicator in 7 Cep 
(the Ca II A8662 emission line index), and thus that the 
velocity variations were spurious and probably due only 
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to changes in the spectral line profiles caused by surface 
inhomogeneities (spots) driven by stellar rotation. More 
recently the planetary interpretation was reinstated by 
Hatzes et al. ( 2003) on the basis of new high-precision ve- 
locity observations at the McDonald Observatory. They 
showed convincingly that the 2.5-yr variation is coherent 
in phase and amplitude throughout the entire 20-yr in- 
terval covered by the merged CFHT and McDonald data 
sets, as would be expected for Keplerian motion, and that 
no changes were observed in the spectral line bisectors. 
On the other hand, a careful re-an alysis of the changes 
in the activity indicator reported bv IWalker et al.l ()1992D 
revealed that the periodicity of the Ca II A8662 measure- 
ments (2.14 yr) is not only slightly different from that in 
the velocities, but it is transitory in nature, thus ruling 
out a connection. 

7 Cep also carries the distinction of being among the 
first planet host stars to be found in a binary system, 
which raises interesting issues related to the dynam- 
ical staHlit^;_of_such^ A recent study 
by iRaghavan et al.l l)2006'l points out that among the 
known planet host stars 7 Cep happens to be the sys- 
tem with the smallest ratio (^11) between the size 
of the binary orbit and the planetary orbit. How- 
ever, the outer orbit is at present poorly known, and 
the secondary star is presumably very faint and has 
never been seen. Reported values for the binary pe- 
riod have ranged between 29.9 yr (tW alker e t alJll99^ 
and 66 yr ( Griffin, Carauillat &: Ginest et 2002), and have 
been based on only part of the data available, in some 
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cases spanning much less than a full cycle. A number of 
authors have carried out numerical investigations of the 
gravitational influen ce of the secondary star on th e orbit 
of the planet (e.g., |Dvorak et al.ll200a fThcbault^TaO 
120041 iHaghighiDOur 2006), but have used rather differ- 
ent parameters for the binary or have pointed out the 
uncertainty in those elements as a limiting factor. 

The motivation for this paper is thus threefold: i) To 
improve the determination of the orbit of the secondary 
star (including for the first time an estimate of the incli- 
nation angle, and of the mass of the secondary) in order 
to allow more definitive dynamical studies of the stabil- 
ity and evolution of the system. We do this by using 
all available radial velocity data for 7 Cep including new 
measurements reported here and other historical obser- 
vations not previously used. We incorporate also astro- 
metric measure ments from the Hipparcos mission ("ab- 
scissa residuals" ; lESAII1997l) as well as transit circle and 
other positional information spanning more than a cen- 
tury, a) To carry out a critical review of previous stud- 
ies of the physical properties of the primary star and use 
all available information to estimate its absolute mass, 
a key parameter influencing the mass of the substellar 
companion, iii) To place firm dynamical upper limits on 
the mass of this companion by taking advantage of the 
high-precision Hipparcos intermediate data and model- 
ing the reflex motion of the primary star on the plane of 
the sky. We show that this modeling allows us to con- 
flrm the substellar nature of the companion, although it 
is not yet possible to rule out a mass in the brown dwarf 
regime. 

2. OBSERVATIONAL MATERIAL 

We describe here all spectroscopic and astrometric 
measurements of 7 Cep of which we are aware that have 
a bearing on the motion of the star, with the goal of 
combining them into a global orbital solution in Sj^l 

2.1. Radial velocities 

The high-precision Doppler measurements of 7 Cep 
have been described in detail by Hatzcs et al. (2003). 
They consist of 4 separate data sets corresponding to 
different instrument conflgurations: three from the Mc- 
Donald Observatory (referred to below as McDonald I, 
McDonald II, and McDonald III, following 'H atzes et al.l 
|2003), and_ one from t he CFHT, which is the data set 
of IWalker et all 1)1992(1 . The nominal precision of these 
measurements ranges from about 8 ms~^ to ~30 ms~^, 
and they are all differential in nature as they rely on the 
use of telluric O2 lines as the velocity metric, or on lines 
of hydrogen fluoride or iodine gas that play the same 
role. Taken together these velocities cover the interval 
1981.4-2002.9 , which includes per iastron passage in the 
binary orbit. iHatzes et al.l l|2003|) combined these data 
and solved simultaneously for the outer orbit and the or- 
bit of the planet. The time span of the observations is 
less than half of their estimated binary period of 57 yr. 

Beginning in the late 1970's 7 Cep was moni- 
tored spectroscopically using m ore traditional means by 
iGriffln. Carauiflat fc Gii^este^ pOOl . To their own ob- 
servations with several different instruments in Cam- 
bridge (England), Haute-Provence (France), and Victo- 
ria (Canada) , they added a subset of the high-precision 
velocities mentioned above as well as other velocities col- 



lected from the literature in an effort to extend the time 
coverage and better constrain the outer orbit. These 
include measures pubhshed bv iBeavers fc Eitteil l|1986D 
made in 1978-1980, and most importantly the veloc- 
ities obtained at the Li ck Observatory in 1896-1921 
ijCamnbell fc Moorell928D. All these mea s ureme nts were 
placed bv IGriffln. Carauillat fc: Ginesteti l|2002() on the 
same zero point (corresponding to their Cambridge in- 
strument), and have formal uncertainties ranging from 
0.2 kms^-'^ to 0.9 kms"-'^. We adopt these 77 measure- 
ments as published. These authors noted an unfortunate 
gap of some 50 years in the velocity coverage for 7 Cep 
that complicates the determination of the orbital period 
(see also ti3.1l) . In order to distinguish between two pos- 
sible periods (66 yr and 77 yr) allowed by the radial ve- 
locity data they used, they considered also other mea- 
surements from the lit e rature in the inte r val 1902-1907 
llFrost fc Adamal IT903t IBelopolskvl IT90I ISlipheJ ITaol 
lKiistneilll908l ^. Thev found that those velocities favored 
the 66-yr period, although they did not actually make 
use of them in their orbital solution because of their un- 
certain zero point. 

Our own contribution to the observational material is 
twofold. On the one hand we have derived 3 new ve- 
locities for 7 Cep based on archival spectra collected at 
the Harvard-Smithsonian Center for Astrophysics (CfA) 
using an echelle spectrograph on the 1.5-m Tillinghast 
reflector at the F. L. Whipple Observatory. The nominal 
precision of these measurements is around 0.3 kms^^ for 
a bright and sharp- lined star such as this. For details 
on the reduction procedures we refer the reader to the 
description by iTorres. Neuhauser fc Guentherl l(200^ . 
While this contribution is modest by comparison to the 
material described earlier, it does extend the time cov- 
erage to the end of 2004, and the velociti es are on 
a wcll-deflned system (s ee Stcfanik. Latham 'fc Torre3 
ri999tlT,atha,m et a,l.ll200l. 

Given the poor spectroscopic coverage prior to 1978, 
we carried out a careful search of the literature for ad- 
ditional measurements that might help constrain the 
outer orbit. Aside from the 1902-1907 sources men- 
tio ned above that were used only a s supp orting evidence 
by iGriffin. Carauillat fc Ginesteij l)2002tl . a number of 
other velocity sources were found but their zero points 
are generally unknown, so the measurements cannot be 
combined at face value. Thus as a second contribution 
we relied on the extensive CfA database of ^250,000 
spectra to place all of those scattered measurements of 
7 Cep onto a uniform frame of reference. This was ac- 
complished by using measurements for other stars also 
reported in each of these sources, and comparing them 
with newly derived velocities for those same "standards" 
from CfA spectra obtained at one time or another over 
the past 25 years. Details of this procedure are pro- 
vided in the Appendix. Of p articular relevance are the 
7 Cep measurements by |K jg;rgaard et al. (1981) made 
injL977. ISm^wden fc Young ( 20051 in 1972-1974, BoiilonI 
(1957^ in 1955, and Harper \l9^ in 1921. The preci- 
sion of those velocities ranges from 0.5 kms^^ to about 
1.9 kms""'^. We list them in Table ^a-long with our own 
measurements, all on the CfA system. 

Despite our attempts to establish their zero points, 
two of the sources of historical velocities showed large 
discrepancies when compared with other data taken at 
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similar times, or pre sented othe r prob lems. The series 
of measurements by iBelonolskvl l|1904l) contains only 5 
other stars usable as standards, and the offset required to 
place those velocities on the CfA system has the largest 
uncertainty ('^1 kms~^). The corrected 7 Cep velocities 
from 1903 are some 3 kms~^ too high. Three velo cities 
meas ured at Mt. Wilson Observatory in 1915-1917 ijAbti 
Il973|) show the largest spread of any data set (4 kms~^). 
The zero point of those measurements is very difficult 
to establish because of the variety of instruments and 
telescopes used, which are not always indicated in the 
original publication. The average of the corrected veloci- 
ties for 7 Cep shows a discrepancy of 5 kms""'^ relative to 
others made within a few years. We have therefore not 
made use of either of these two data sets in our orbital 
solution described in O Several high-dispersio n plates of 
7 Cep were obtained bv lKoelbloed fc van Parad iis (1975) 
in 1963-1964, a critical time in the observational his- 
tory of this object, but unfortunately the authors ap- 
pear not to ha ve measured radial velocities. Finally, 
iRucihski &: Stan iucha (.1981) published a velocity mea- 
surement for 7 Cep made in 1979, but all of the other 
stars reported in their paper happen to be variable, and 
so cannot be used as standards. 

2.2. Astrometry 

7 Cep was observe d by the H ipparcos satellite be- 
tween 1989 and 1993 l)ESAII1997|) . These accurate one- 
dimensional astrometric measurements were used by the 
Science Team to derive the position, proper motion, and 
trigonometric parallax of the object (tthip = 72.50 ± 
0.52 mas) as reported in the main catalog. The as- 
trometric solution revealed a measurable acceleration on 
the plane of the sky (proper motion derivatives) in the 
amount of d^a/dt — -1-1.51 ± 1.12 mas yr~^ in Right 
Ascension and a more significant dfis/dt — -1-6.10 ± 
1.11 mas yr^^ in Declination. This acceleration is of 
course due to the binary nature of the object, and was 
accounted for in deriving the parallax. 

As we demonstrate below, the binary motion at the 
epoch of the Hipparcos observations is such that we ex- 
pect some curvature on the plane of the sky that should 
be detectable in the measurements. We have therefore 
made use of these observations (available in the form of 
"abscissa residuals" ) in our orbital solution described be- 
low, since they are complementary to the spectroscopic 
observations and provide new information. A total of 
76 such measurements were obtai ned by the two inde- 
pendent data reduction consortia l)ESAII1997j) . and the 
median error for a single measurement is 1.9 mas. 

Because of the relatively short time span of these 
observations compared to the binary orbital period, it is 
almost certain that part of the orbital motion has been 
absorbed into the proper motion components reported 
by Hipparcos. This is in fact a way in which many 
long-period binaries have been discovered in the past, 
on the basis of the apparent variability of their proper 
motions whe n comp ut ed at different epochs (see, e.g., 
iWielen et aD 119991: iContcharov. Andronova & Titovl 
120001: iMakarov fc Kaolaiil l2005|) . Precise l y this effect 
was pointed out for 7 Cep by iHeintd l|199(1f) . who 
noticed a significant change mostly in fis over several 
decades. Therefore, to make proper use of the Hipparcos 
intermediate data to extract information on the binary 



orbit it is necessary to constrain the proper motion 
by other means in order to model the orbital motion 
without risking systematic errors. 

Initially we considered using the proper motion fo r 
7 Cep reported in the Tycho-2 catalog fH0s: et al. 2000a) , 
which relies on ground-based positional measurements 
made over many decades and is constrained at the re- 
cent epoch by the Tycho-2 position. This long baseline 
presumably averages out any perturbations due to or- 
bital motion if the period is significantly shorter than 
this. The Tycho-2 proper motion is in fact quite dif- 
ferent from the Hipparcos determination, which is effec- 
tively "instantaneous" at the mean epoch ~1991.25. In 
the case of 7 Cep, however, the orbital period is not negli- 
gible compared to the time span, and we were concerned 
that fia and might be biased. Evidence that the or- 
bital motion is detectable in the individual positional 
measurements from tra nsit circle observations wa s indeed 
presented by IContcha rov. Andronova fc Titovl (|2000), 
who inferred from them a period of about 45 yr for 
the binary. We therefore chose to make use of the 
individual positions from ground-based catalogs going 
back to 1898, kindly provided by S. Urban of the U.S. 
Naval Observatory (USNO). Additional measurements 
from the Carlsberg Meridian Catalogs (CMC; see, e.g., 
ICar lsberg Meridian Catalogue 1989) were provided by 
G. Gontcharov (Pulkovo Observatory) or obtained from 
the literature. All of these measurements have been 
reduced to the International Celestial Reference Frame 
(ICRF), effectively represented in the optical by the Hip- 
parcos catalog, and their nomin al precision v aries be- 
tween about 50 and 500 mas fsee lHog et aHMOQlh . We 
Hst them in Table [3 

3. ORBITAL SOLUTION 

The combination of the radial velocity mesurements 
and the astrometry makes it possible to derive the com- 
plete set of elements describing the binary orbit in 7 Cep. 
The inclination angle is of particular interest because 
when combined with the spectroscopic mass function it 
provides the information needed to compute the mass 
of the secondary star, given an estimate of the primary 
mass. The substellar companion to the primary intro- 
duces additional components of motion that we model 
simultaneously. Given that the outer orbit is an order of 
magnitude larger than the inner orbit (see to first 
order we assume here that they are decoupled, i.e., that 
the outer one may be treated as corresponding to a "bi- 
nary" composed of the secondary star (B) and the center 
of mass of the inner pair (A). Orbital elements that refer 
to the outer orbit are indicated below with the subindex 
"AB" , and those pertaining to the inner orbit are dis- 
tinguished with a subindex "A" . The primary star itself 
is referred to as "Aa" following the traditional spectro- 
scopic notation, and the planet (indistinctly called also 
'substellar companion') as "p" , for simplicity. 

The radial velocities allow us to solve for the period, 
center-of-mass velocity of the (triple) system, eccentric- 
ity, velocity semi-amplitude, longitude of periastron, and 
time of periastron passage in the outer orbit: {Pab, 
7, CAB, Ka, CO a, Tab}- The high-precision velocities 
constrain the spectroscopic elements of the inner (plan- 
etary) orbit: {Pa, ga, Ka^, tfAa, Ta}- Because the 
high-precision velocities are differential, an offset must 
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be determined to place them on the frame of the abso- 
lute velocities, for which we have chosen the Griffin data 
set as the reference. We therefore solved for 4 additional 
parameters representing these offsets, one for each data 
set: ARVi, Ai?V2, and Ai?V3 for the groups referred 
to as McDonald I, McDonald II, and McDonald III (see 
323), and ARVi for the CFHT data set. The CfA ve- 
locities and other historical data sets placed on the CfA 
system were considered as a single group, and one ad- 
ditional parameter ARV5 was included to represent the 
shift relative to Griffin. 

Preliminary estimates suggested the secondary star is 
very small compared to the primary, and we may assume 
here that it contributes no light. The Hipparcos obser- 
vations therefore refer strictly to the primary as opposed 
to the center of light, and provide a constraint on the 
orientation of the outer orbit (inclination angle «ab and 
position angle of the ascending node ^Iab, referred to 
the equinox of J2000) as well as on the angular scale of 
the orbit of the inner binary relative to the barycenter 
(a^). We show below that the astrometric measurements 
do not, however, resolve the wobble of the primary star 
caused by the planet. We point out also that there are 
no available measurements of the relative position be- 
tween the two stars, since the secondary has never been 
resolved. The use of the Hipparcos measurements in the 
global solution introduces several other parameters that 
must be solved for, including corrections to the catalog 
values of the position of the barycenter ( Aa* , AS) at the 
mean reference epoch of 1991.25, and corrections to the 
proper motion components (A/i* , A^g)^. In principle we 
also need to solve for a correction to the Hipparcos par- 
allax. However, the fact that the spectroscopic elements 
of the outer orbit are solved for at the same time intro- 
duces a redundancy, and the parallax (which in this case 
would be termed an "orbital" parallax) can be expressed 
in terms of other elements as 



Pab Ka\/1- Cab ' 

where the period is given in days and Ka in kms~^. We 
have therefore chosen to eliminate the parallax correction 
as an adjustable parameter. 

By combining complementary observations of different 
kinds the global solution is strengthened. The ground- 
based positional measurements provide the tightest con- 
straint on the proper motion and position of the barycen- 
ter. This breaks the strong correlation between proper 
motion and orbital motion in the Hipparcos observations, 
and enables those measurements to provide a constraint 
on the angular scale and orientation of the outer orbit, 
even though their coverage is only a small fraction of 
the orbital period. Some information on the scale and 
orientation, as well as on the outer period, is provided 
also by the positional measurements, while the veloci- 
ties contribute most of the weight to the period, shape, 
and linear scale of the binary orbit. The elements of the 
inner orbit are constrained only by the high-precision 
velocity measurements, and are only weakly dependent 
on the outer orbit. Light-travel effects in the inner or- 
bit are negligibly small. The formalism for incorporat- 

^ Following the practice in the Hipparcos catalog wo define 
Aa* = Aa cos 5 and A/i* = Afia cos S. 



ing the abscissa residuals from Hipparcos into the fit 
follows clos ely that described by fvan Leeiiwen fc EvansI 
l)1998|) and iPourbaix fc Jorissenl l)200(]|) . including the 
correlations between measurement s from the two inde- 
pendent data reduction consortia fESA'"1997). In using 
the ground-based catalog positions the parallactic mo- 
tion was accounted for in our model, given that the pre- 
cision of some of the more recent measurements is com- 
parable to the parallax. 

Altogether there are 23 unknowns that we solved for 
simultaneo usly, using standa rd non-linear least-squares 
techniques l|Press et al.lll992l p. 650). The solution con- 
verged quickly from initial values of the elements cho- 
sen from preliminary fits or by an extensive grid search, 
and experiments in which we varied the initial values 
within reason yielded the same results. A total of 446 
individual observations were used from 10 different data 
sets, as follows: 107 classical radial velocities (77 from 
Griffin, 30 from CfA and other literature sources), 199 
high-precision velocities (68 from CFHT, 43 from Mc- 
Donald I, 49 from McDonald II, and 39 from McDon- 
ald HI), 76 one-dimensional Hipparcos measurements, 
and 64 ground-based catalog coordinates (split into two 
data sets of 15 and 17 pairs of Right Ascension and Dec- 
lination measurements). Weights were assigned to the 
measurements according to their individual errors. Since 
internal errors are not always realistic, we adjusted them 
by applying a scale factor in such a way as to achieve a 
reduced value near unity separately for each data set. 
This was done by iterations. These scale factors were all 
close to unity for most of the velocity sets, and somewhat 
larger for some of the ground-based catalog positions^. 

The results are given in Table |2| along with derived 
quantities such as the position of the barycenter at the 
mean epoch of the Hipparcos catalog (1991.25), the par- 
allax and proper motion components, and the mass func- 
tion of the stellar binary. Other derived quantities are 
described below. The elements of the planetary orbit 
are not sign i ficant ly different from those reported by 
iHatzes et al.l ()2003t) since that orbit depends essentially 
only on the high-precision velocities, for which we used 
the same data they used. The parallax is also not ap- 
preciably different from the Hipparcos value, although 
our uncertainty is somewhat smaller. The proper mo- 
tion components, on the other hand, are considerably 
different from their catalog values, as anticipated above. 

The radial velocity measurements are represented 
graphically in Figure ^ The top panel shows only the 
classical measurements. The small undulations in the 
computed curve are produced by the wobble of the pri- 
mary star with the 2.47-yr period of the substellar com- 
panion. Although the phase coverage in the outer orbit 
is incomplete, the period of the binary is fairly well es- 
tablished thanks in part to the very high precision of 
the McDonald and CFHT data (see also aTT|l . These 
are shown separately in the lower panel, where the error 
bars are smaller than the size of the points. The reflex 
motion of the primary due to the substellar companion 
is shown as a function of orbital phase in Figurc|21 where 

The scale factors derived are 0.96 (Griffin velocities), 0.93 (Mc- 
Donald I), 1.10 (McDonald II), 1.02 (McDonald III), 1.45 (CFHT), 
0.94 (CfA), 1.62 (USNO Right Ascensions), 1.50 (USNO Declina- 
tions), 0.83 (CMC Right Ascensions), 0.87 (CMC Declinations), 
and 0.86 (Hipparcos). 
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the motion in the outer orbit has been subtracted from 
the individual data sets. The residuals from the new 
CfA velocities and from other values from the literature 
sources placed on the same system are given in Table ^ 

The path of 7 Cep on the plane of the sky is repre- 
sented in Figure |3| where the axes are parallel to the 
Right Ascension and Declination directions. The solid 
curve is the result of the contributions from the annual 
proper motion (arrow), the parallactic motion, and the 
motion in the binary. The wobble due to the substellar 
companion is negligible on the scale of this figure (see 
also The predicted location of the one-dimensional 
Hipparcos observations is indicated by the dots on the 
curve. As stated earlier, their typical uncertainty is 1.9 
mas. For illustration purposes, the dotted line in the 
figure starting at the location of the first Hipparcos ob- 
servation shows the path the star would follow without 
the perturbation from the orbital motion in the binary. 

The orbit of the primary star in 7 Cep around the cen- 
ter of mass of the binary is shown in Figure 01 with a 
semimajor axis of about 325 mas. The direction of mo- 
tion is retrograde (arrow). The intersection between the 
orbital plane and the plane of the sky (line of nodes) is 
represented by the dotted line. The section of the orbit 
covered by the Hipparcos mission is indicated with filled 
circles, and the open circle labeled "P" represents peri- 
astron. A close-up of the area around the Hipparcos ob- 
servations is shown in Figure |S1 Because these measure- 
ments are one-dimensional in nature, their exact location 
on the plane of the sky cannot be shown graphically. The 
filled circles represent the predicted location on the com- 
puted orbit. The dotted lines connecting to each filled 
circle indicate the scanning direction of the Hipparcos 
satellite for each measurement, and show which side of 
the orbit the residual is on. The short line segments 
at the end of and perpendicular to the dotted lines in- 
dicate the direction along which the actual observation 
lies, although the precise location is undetermined. Oc- 
casionally more than one measurement was taken along 
the same scanning direction, in which case two or more 
short line segments appear on the same dotted lines. 

The motion of 7 Cep A in the binary orbit is discernible 
in the ground-based catalog measurements taken over the 
last century, although only in the Declination direction. 
This is illustrated in Figure El The amplitude of motion 
in the R.A. direction is much smaller because of the ori- 
entation of the orbit, which is mostly North-South. The 
more recent measurements since 1980 are much more pre- 
cise. That section of the orbit is shown on a larger scale 
in Figure The residuals of all ground-based measure- 
ments from our orbital solution are given in Table (21 

3.1. The constraint on the binary period 

The poor observational coverage of 7 Cep prior 
to 1980 has made it difficult to establish the pe- 
riod of the outer orbit in previous studies, particu- 
larly since th e secondary star has never been resolved. 
I Walker et al.l |l992,) gave a rough estimate of 29.9 yr 
based on only 10 yr of high-precision veloc ity cover- 
age. iGontcharov. Andronova fc Titovl l)2000() used the 
ground-based catalog positions spanning a little less 
than six decades and inferred a period of 45 yr. A 
re-reduction of those same data (G. Gontcharov 2006, 
private communication) does not show that periodic- 



ity as clearly, however. The radia l- velocity study by 
IGriffin. Carauillat fc Ginestetl lj2002|) took advantage of 
some of the historical measurements going back more 
than 100 yr, and found that a 50-yr gap in the data near 
the middle of the last century allowed two possible pe- 
riods giving fits of similar quality: '^66 yr and ~77 yr. 
On the basis of other ob servational evidence they chose 
the short orbital period. iHatzes et al.l l)2003|) used 20 yr 
worth of high-precision velocity measurements and de- 
rived a period of 57 yr. 

The simultaneous use of all of the above measurements, 
and the addition of other observations (including more 
recent velocities as well as historical velocities, and the 
Hipparcos measurements), have allowed us to finally con- 
strain the binary period without ambi guity to a value 
of 66.8 yr, thus proving iGriffin. Carauillat fc Ginestetl 
(2002^) essentially correct. To illustrate the improvement 
brought about by the added observations, we have re cre- 
ated the fit by IGriffin. Carauillat fc Ginestetl l|2002(l by 
using the same set of observations they used'^, ignoring 
the velocity perturbation from the substellar companion, 
as they did. In FigurejSlwe show the reduced of the fit 
for a range of fixed orbital periods, with the remaining 
orbital elements adjusted as usual to minimize x^. We 
then repeated this exercise using the data that went into 
our own solution, this time accounting properly for the 
planetary compani on. The dashed curve corresponding 
to the solution bv IGriffin. Carauillat fc Ginestetl l)2002f) 
shows two local minima at ~66 yr and ~77 yr, as found 
by those authors. The solid curve corresponding to the 
solution in this paper that includes all available obser- 
vations has a single minimum at 66.8 yr. The formal 
uncertainty in this value is 1.4 yr, or 2%. 

4. PHYSICAL PROPERTIES OF THE PRIMARY AND 
SECONDARY STARS 

In order to take full advantage of the orbital solution 
presented above and estimate the mass of the unseen sec- 
ondary star, as well as to place limits on the mass of the 
substellar companion, we require an estimate of the mass 
of the primary star itself (AfAa). Mass estimates in the 
literature for 7 Cep have varied by more than a factor 
of two (between ^^0.8 Mq and ~1.7 Mq), which is some- 
what surprising for such a bright and well-studied star 
but may perhaps be explained by its present evolution- 
ary state and other uncertainties (see below). We wish 
to constrain it to much better than this to avoid prop- 
agating the uncertainty to other quantities that depend 
on the mass. In this section we therefore examine the 
available observational material carefully and critically, 
making use of current stellar evolution models to arrive 
at the best possible estimate for Mab.- We discuss some 
of the other estimates as well in an attempt to under- 
stand the differences. 

The brightness of 7 Cep has made it an easy target for 
spectroscopic studies to determine both the effective tem- 
perature and chemical composition of the star. These, 
along with other properties, are essential in order to es- 
timate its absolute mass. In Table 01 we have collected 

^ The only difference in our recreated solution is that we used 
the high-precision velocity measurements as published, whereas 
Griffin. Carauillat & Ginestet (2002) used values read off from a 
figure, since some of the measurements had not yet been reported 
in tabular form in the literature. 
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the results of nearly two dozen separate investigations 
carried out over the past 40 years. We consider here 
only determinations of Toff and [Fe/H] that are purely 
spectroscopic. Our own temperature estimate from the 
spectra described in i )2.1l is listed as well (see [0for the 
details of our procedures). For the most part the 22 
independent metallicity determinations show reasonable 
agreement within the errors, and yield a weighted average 
of [Fe/H] ~ +0.01 ± 0.02, or very nearly solar. Further 
comments on this value are given below. The weighted 
average effective temperature is Tcs = 4852 ±26 K from 9 
spectroscopic measurements including our own. The un- 
certainties given here are statistical errors that account 
for the different weights as well as the scatter of the indi- 
vidual [Fe/H] and Tcff measurements, but not for possible 
systematics. In the following we adopt for these averages 
more conservative errors of 0.05 dex and 100 K, respec- 
tively. 

Temperature estimates for the star have been derived 
on numerous occasions also from color indices in a variety 
of photometric systems. In order to bring homogeneity 
to this information we have compiled the available pho- 
tometry for 7 Cep in eight different systems (Johnson, 
Stromgren, Vilnius, Geneva, Cousins, DDO, 2MASS, 
and Tycho) mostly from the photometric database main- 
tained by Mermilliod. Mcrmilliod & Hauck (1997), and 
we have used t he color / temperature c alibra tions for gi- 
ant stars from iRamfrez fc MelendezI Emii for 13 dif- 
ferent photometric indices. Interstellar reddening has 
been ignored here in view of the close distance to the 
star (13.8 pc), but we have accounted for the very small 
metallicity correction in each calibration based on the 
discussion above. The results are collected in Table 
where the uncertainty of each temperature estimate in- 
cludes the contribution from photometric errors as well 
as the statistical uncertainty of the calibration, added in 
quadrature. The weighted average of these determina- 
tions is Toff — 4754± 17 K, although we prefer 100 K as a 
more realistic error to account for unquantified system- 
atics. The spectroscopic and photometric temperature 
estimates thus differ by only 100 K, and we adopt here 
the compromise value of Tcff — 4800 ± 100 K. 

Two additional properties of the star that can 
be determined very accurately are the absolute vi- 
sual magnitude and the linear radius. The abso- 
lute mag nitude fo llows from V = 3.213 ± 0.007 
ijMermilliod. Mermi lliod fc Haucldll997^ and our paral- 
lax for the system (TableEJ, and is My = 2.521 ± 0.014. 
The angular diameter of 7 Cep has been measured di- 
rectly with high precision bv iNordgren e t eJ] |)1999D us- 
ing the Navy Prototype Optical Interferometer, and is 
(j) = 3.24 ± 0.03 mas (limb-darkened value). Combined 
once again with our parallax, this measurement yields 
the linear radius a,s R — 4.790 ± 0.052 R©, which has a 
formal precision just over 1%. 

In Figure |51 we compare the measured temperature 
and absolute visual magnitude of the star with e volu- 
tionary tracks fr om the series by iiYi et all IpOOl and 
iDemaraue et all (|2p04) for the composition established 
above. Tracks are labeled with the mass in solar units. 
The star is seen to be in the first ascent of the giant 
branch. The shaded error box is shown more clearly 
in the inset, which suggests a mass for 7 Cep slightly 
over 1.2 M© and an uncertainty in that value determined 



almost entirely by the temperature error at this fixed 
metallicity. If the radius is used instead of the tempera- 
ture, the constraint on the mass is considerably improved 
because of the smaller relative error of R. This is shown 
in Figure irm which indicates a mass also close to 1.2 M©, 
consistent with the previous figure. In both cases the un- 
certainty in the metallicity is also important, as a change 
in [Fe/H] shifts the tracks essentially horizontally. 

The optimal value of the mass is one that yields the 
best simultaneous match to the four measured quanti- 
ties (Toff, [Fe/H], Mv, and R) within their stated errors. 
To determine this value, as well as its uncertainty, we 
computed by interpolation evolutionary tracks in a fine 
grid for a range of masses and also a range of metallic- 
ities within the observational uncertainty of [Fe/H]. At 
each point along the tracks we compared the predicted 
stellar properties with the measurements, and recorded 
all models that agree with the observations within their 
errors. All such models are displayed in Figure ^3 in 
a mass/age diagram. It is seen that at each mass the 
range of allowed ages is very narrow. The best match 
is for a mass of Mah — 1.18to ii M©, and the corre- 
sponding evolutionary age is 6.6'!iQ'^ Gyr. All four mea- 
sured quantities are reproduced to well within their er- 
rors (better than 0.3(t), an indication that they are mu- 
tually consistent. The surface gravity predicted by the 
best model is log g — 3.15. An independent age esti- 
mate was obtained by iSaffe. Gomez fc Chavero| l|2005f) 
based on the chromospheric activity indicator log R'n^^ = 
—5.32. Their result (6 . 39 Gy r) using the calibration by 
iRocha-Pinto fc Maciell l)1998|) agrees very well with ours 
formally, although chromospheric ages for older objects 
tend to be rather uncertain^. 

There are significant differences between our mass and 
other recent estimates. Almost all of them rely on evo- 
lutionary models and use different co mbinations of ob- 
servational constraints. For example, iFuhrmarml l|2004fl 
derived a value of 1.59 M© with a formal error less than 
10%, from a fit to the effective temperature and bolomet- 
ric magnitude (derived using the Hipparcos parallax) for 
a fixed metallicity that is higher tha n ours (see Table IH) . 
This mass estimate was adopted bv lHatzes et all 1)20031) 
to infer the minimum mass of the substellar companion 
to 7 Cep. pVffcr et al. (2005) obtained an even larger pri- 
mary mass of 1.7 M© (no uncertainty given) from a fit 
to their own Toff and My (also based on the Hipparcos 
parallax) , using t heir [Fe/H] determination that is again 
higher than ours. lAllende Prieto fc Lamberil l|1999t) used 
My and B ~ V directly and inferred Mas, = 1.33 M©, 
but apparently made no use of any measured metallic- 
ity. A lowe r mass than ours (A/a a, — 1.0 ± 0.2 M©) was 
derived by iLuck fc Challener (1995) from the luminos- 
ity and temperature they determined for 7 Cep, along. 
with t heir [Fe/H] value, which is close to solar. Glcbocy 
l)1972l) obtained 1.5 M© employing a similar method, but 
adopted a metallicity much lower than ours. Except for 
the latter study, the evolutionary models used by most 
of these authors are similar enough that the differences 
in mass must be due in large part to the observational 

* An additional age estimate bv|[yaffe, G omez Chaverd J200F)I1 
based on the calibration by iDonahug^^olil) gave the value 
14.78 Gyr for 7 Cep, which, however, is older than the age of 
the Universe. 
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constraints, particularly the temperature and metallicity. 
We note also that none of these studies have made use of 
the measured angular diameter of the star, which appears 
to be very accurate. An entirely different approach was 
followed by iGratton et al.l l)1982l) , who inferred a mass 
from their spectroscopic log 17 determination along with 
a linear radius derived from surface brightness relations 
[R = 6.5 R0). Their value is MAa = 0.89 ± 0.19 Mq. 

The larger estimates of tend to be those based on 
hotter temperatures and also higher metallicities. Fur- 
thermore, a look at Table0]shows that while most of the 
metallicity determinations for 7 Cep are close to solar, 
all the higher values have been reported only in the last 
few years, and they tend to go together with hotter tem- 
peratures (see Figure I12II . We note in this connection 
that 7 Cep has been considered a member of a group of 
evolved stars displaying CN bands that are stronger than 
usual ('strong- CN stars ' or 'v e ry-strong-lined stars'; see, 
e.g., J Spinrad fc Tavloil 119691: iKeenan. Yorka fc WilsonI 
It was classified by the latter authors as Kl III- 
IV CN 1, which corresponds only to a marginal strong- 
CN star. These objects have had a controversial history, 
occasionally having been considered to be super-metal- 
rich. Other studies have disputed this, however. For ex- 
ample, some of the giants in the open cluster M67 have 
been found to have strong CN features even though the 
chemical composition o f this cluster is believed to be es- 
sentially solar (see, e.g.. lLuck fc Challeneilll995|) . We re- 
fer the reader to the latter work (and references therein) 
for an excellent summary of the subject and a list of 
possible explanations for the CN phenomenon. 

As tempting as it may be to place higher confidence in 
some of the more recent [Fe/H] studies that have found a 
metal-rich composition for 7 Cep, it is difficult to ignore 
the large body of equally careful determinations yielding 
a comp osition closer to solar. This includes the recent 
work of iFranchini et all l|2004f) . which not only gives a 
slightly subsolar metallicity but also happens to have 
the smallest formal uncertainty; their result is [Fe/H] 
= —0.066 ± 0.034. As a test we repeated the com- 
parison with stellar evolution models described earlier, 
but adopting the spectroscopic temperature and metal- 
licity determinations of each of the recent studies that 
give super-solar abundances. In no case did we find a 
model that is simultaneously consistent with all four of 
the quantities Tcff, [Fe/H], My and R, within their un- 
certainties. We are led to conclude, therefore, that the 
chemical composition of 7 Cep is not significantly higher 
than solar, and we adopt in the following the mass we de- 
termined above with the most conservative of the asym- 
metric error bars: AfAa — 1-18 ± 0.11 M0. 

With this value and the mass function from our or- 
bital solution the mass of the unseen stellar companion 
is Mb = 0.362 ± 0.022 M©, where the error is computed 
from the full covariance matrix resulting from our fit (in- 
cluding cross-terms) and accounts also for the primary 
mass uncertainty, which represents the dominant contri- 
bution. Thus the secondary is most likely a late-type 
star^ of spectral type approximately M4. The angular 
semimajor axis of the relative orbit between the primary 

^ For completeness we mention here two alternate possibilities, 
although we consider them much less likely. One is that the sec- 
ondary is a white dwarf. In this case its low mass would make it 
a helium-core white dwarf, which are the products of binary evo- 



and secondary becomes a^g — 1. 382±0. 047 arcsec, which 
corresponds to 19.02 ± 0.64 AU. 

With the secondary mass known, it is of interest to 
compute its brightness relative to the primary in or- 
der to assess the chance of detecting it directly, most 
likely in the infrared. The brightness measurements of 
the primary itself in the near infrared are rather un- 
certain because the star saturated the 2MASS detec- 
tors (see Table EJ. From our best model fits we derive 
absolute magnitudes of Mh — 0.28 and Mk — 0.19 
in the Johnson system, which are actually consistent 
with the values inferred from the 2MASS photometry 
within their large errors. The brightness of the sec- 
ondary star may be estimated also from stellar evolu- 
tion models. For this we have used the calc u lation s by 
Baraffe et alJ (|!99RV since those of iYi et all l|200lD are 
not intended for low-mass stars. For the age we estab- 
hshed above we obtain My — 10.92, Mh ~ 6.83, and 
Mk ~ 6.56, whi ch we have plac ed on the same photo- 
metric system as 'Yi et al. (2001) following the prescrip- 
tion by Bessell fc Brett (1988). Thus, the secondary is 
expected to be ~8.4 mag fainter than 7 Cep A in V, ~6.6 
mag fainter in H, and '~6.4 mag fainter in K.^ 

The orbital elements in Table El allow the relative po- 
sition of the unseen secondary to be predicted. We 
note, however, that the scale of the relative orbit still 
depends critically on the assumed primary mass as 
(MAa + Me)^/^, in which the secondary mass itself scales 
as (AfAa + Mb)^^'^. As seen earlier MAa is quite sensitive 
to the adopted temperature and metallicity. A dynami- 
cal (hypothesis-free) estimate of the masses of both stars 
and a direct measure of the semimajor axis a^g will be 
possible once 7 Cep B is detected and its path around the 
primary measured over at least a portion of the orbital 
cycle. 

5. THE MASS OF THE PLANETARY COMPANION 

The reficx motion of the primary star along the line 
of sight in response to the putative substellar companion 
leads to a mass function of f{Mp) = (1.83 ± 0.32) x 
10~^ Mq from our orbital fit. With the adopted value of 
MAa this corresponds to Mp sin za — 1.43 ± 0.13 Mjup, 
which is only slightly smaller than the valu e Mn s inzA = 
1.7 ± 0.4 Mjup reported bv iHatzes erahl l|2003(l . The 
difference is due almost entirely to the choice of primary 
mass, for which they used MAa = 1-59 Mq. 

The perturbation on the primary star on the plane 
of the sky caused by the substellar companion is ex- 
pected to be small, although it depends obviously on the 
unknown inclination angle «a (and through it on Mp). 
Given that the Hipparcos measurements are fairly pre- 
cise, we attempted to determine this astrometric wobble 

lution involving mass transfer through Roche-lobe overflow. Not 
only is it difficult to see how the substellar companion could have 
survived in this environment (unless it formed later, perhaps from 
remnant material), but there also appears to be no evidence of a 
(presumably hot) white dwarf in ultraviolet spectra of 7 Cep. The 
other possibility is that the companion is itself a closer binary com- 
posed of smaller main-sequence stars. In this case their combined 
brightness would be significantly less than that of a single M4 star 
of the same mass, making it more difficult to detect 7 Cep B. 

^ The brightness of the secondary in V may be overestimated 
by up to 0.5 ma g due to the possibility of missing opaci ties in the 
models (see, e.g.. lDelfosse et al.l200(l : )Chabrier et aLKOOfiTl . which 
would affect the optical the most. 
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simultaneously with the other elements by incorporating 
additional adjustable parameters into the model. Four of 
the elements of this astrometric orbit are already known 
from spectroscopy (Pa, ca, WAa, and Ta)- The remaining 
three are the angular scale (semimajor axis) of the orbit 
of the primary around its center of mass with the planet 
("^Aa)' inclination angle of the planetary orbit (ja), 
and the position angle of the ascending node (f^A, J2000). 
Since spectroscopy gives the projected linear semimajor 
axis (oAa sin ia), and the parallax is a known function of 
other elements (see eq. we take advantage of the re- 
dundancy to eliminate the angular semimajor axis a'J^^ as 
an adjustable parameter, given that it can be expressed 
as 



^Aa — '^A 



Pa Ka, 
Pab Ka 



- e\ smtAB 
CAB siniA 



(2) 



A solution with a total of 25 adjustable parameters did 
not yield a statistically significant detection of the astro- 
metric wobble: the best fit corresponded to an inclina- 
tion angle 19° from face-on, implying a semimajor axis 
OAa = 0-46 ± 0.36 mas and a planet mass around 4 Mjup. 

In order to place a meaningful upper limit on iWp we 
explored the full range of possible values of ?a and 17a to 
identify the area of parameter space where the solutions 
become inconsistent with the observational errors. For 
each pair of fixed values of ia and VLa we solved for the 
other 23 parameters of the fit as usual. A false alarm 
probability can be attached to the A^^ (increase in 
compared to the minimum) associated with each of these 
solutions. In this way we may determine the minimum 
value of siniA (highest value of Mp) for a given confi- 
dence level. This is illustrated in Figure 1131 where we 
show the region of parameter space in the two variables 
of interest along with confidence contours. The light gray 
area corresponds to solutions that can only be ruled out 
at confidence levels up to Icr (~68%), and includes our 
best fit mentioned above (indicated with a plus sign). 
The middle gray area is the region between Icr and 2ct, 
and the dark gray area corresponds to confidence levels 
between 2a and 3cr. At the 2cr level (^^95% confidence) 
the observations rule out companion masses larger than 
13.3 Mjup (or inclination angles less than 6?2 from face- 
on), which would induce reflex motions on the primary 
with a semiamplitude of at least 1.5 mas. This mass cor- 
responds roughly to the conventional boundary between 
planetary and brown-dwarf masses. At a higher confi- 
dence level of 3(7 (99.73%) the mass limit is 16.9 Mjup 
(or lA > 4?9), which would produce a wobble with a 
semiamplitude of about 1.8 mas. There is little doubt, 
therefore, that the companion is substellar. 

6. DISCUSSION AND CONCLUDING REMARKS 

In their paper iHatzes et all l)2003|) attempted to place 
limits on the mass of the substellar companion in a dif- 
ferent way by using their measured projected rotational 
velocity for 7 Cep (wsini = 1.5 ± 1.0 kms~^) along 
with the period they determined for the variation of the 
Ca II A8662 emission-line index (781 ±116 days) and 
the estimated radius of the star (i? = 4.66 R©, adopted 
from iFuhrmann 2004). They relied on two assumptions: 
that the spin axis of the star is parallel to the axis of 
the planetary orbit, and that the period of variation of 
the Ca II index represents the true rotation period of 



the star. The comparison between the measured usini 
and the expected equatorial rotational velocity i^v^a) then 
gives limits on sin i (or sin ?a in our notation) . As it turns 
out, however, there is a mathematical error in their cal- 
culation of Vqo,: they reported iigq = 4.9 kms~^, while 
the correct value is 0.3 kms~^ (see also I Walker et alJ 
Since this is smaller than their vsini, no lim- 
its can be placed on siniA in this way. Their statement 
on the probable mass range of the planetary compan- 
ion is therefore not valid. Either the measured wsini 
is overestimated, or the period of the Ca II variations 
is not the true rotation period of the star. The former 
explanation is perhaps supported by a measurement by 
^rav & Nagar ( 1985), who gave v sini = 0.0±0.8 kms~^ 
(along with a sizeable radial-tangential macroturbulence 
of Crt = 4.2 ±0.6 kms~^, which can effect rotational ve- 
locity m easurements if not properly ac counted for). The 
study bv lde Medeiros fc Mavoil l)1999|) reported v sini < 
1.0 kms~^. Alternatively, the rotation period would 
have to be considerably shorter than 781 days («100- 
500 days), and another explanation would have to be 
found for the variations in the emission-line index. Our 
dynamical constraint on Mp thus shows for the first time 
that the companion is substellar in nature, although a 
mass in the brown dwarf regime (as opposed to the plan- 
etary regime) cannot be completely ruled out with the 
present observations. 

7 Cep is one of more than two dozen examples of sub- 
stellar companions f ound in stellar binaries (see, e.g., 
iRagh avan et al . 2006) . Such systems have attracted con- 
siderable interest in recent years, and numerical studies 
have been carried out specifically for the case of 7 Cep 
to assess not only the dynamical s tability of the orbit 
of the substellar c ompa i iion (e.g., iDvorak et al.l l2003t 
iSolovava fc PittichI I200I lHaghighipo ur"2006'). but also 
the stability of the orbits of other (possibly Earth-like) 
planets that mig ht be present in the h abitable zone of 
the primary star. iThebaiilt et"all l|2004() have also inves- 
tigated the conditions under which the substellar com- 
panion may form through core accretion in the binary 
environment. With a relative semimajor axis for the 
planet orbit of aAa-p = 

1.94 ± 0.06 AU (for an adopted 
primary mass MAa = 1-18 ± 0.11 Mq), the size of that 
orbit is only 9.8 times smaller than the size of the binary 
orbit (19.02 AU; see TableOl), currently the lowest value 
among the known exoplanets in binaries^. Orbit stability 
depends quite strongly on the parameters of the binary 
system, in particular the semimajor axis and eccentricity, 
as well as on the masses of the components. The dynam- 
ical studies mentioned above have all had to make do 
with the rather poorly determined binary properties and 
a lso often inconsisten t results from various authors. 

iHolman fc Wieger^ l)1999() have derived a simple em- 
pirical formula for computing the maximum value of the 
semimajor axis of a stable planetary orbit ("critical" 
semimaj or axis, arnt) in a cop lanar S-type planet-binary 
system. iHaghighipouil lf2006) pointed out in his study 
that the uncertainty in the binary orbital elements made 
for a very large parameter space to be explored numer- 

^ The slig htly smaller orbit size ratio compared to the value of 
~11 given bv lR.aghavan et al.l 12003^ is largely due to the significant 
improvement in the elements of the binary orbit in the present 
work. 
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ically for 7 Ccp. Furthermore, the inclination of the bi- 
nary orbit was unknown at the time and therefore so was 
the mass of the secondary star. As a result, he was only 
able to provide a rather wide range of critical semima- 
jor axes as a function of the adopted binary eccentricity 
(see his Figure 1). With the present study that situ- 
ation has changed, and the critical semimajor axis can 
now be computed directly with a relatively small formal 
uncertainty. We obtain acrit — 3.61 ± 0.36 AU, which is 
considerably larger than the semimajor ax;is of the planet 
orbit, implying the latter is stable if coplanar with the 
binary. 

The combination of classical as well as high-precision 
radial velocity measurements of 7 Ccp with ground- and 
space-based astrometry have allowed a significant im- 
provement in the binary orbital elements (and a first 
determination of the inclination angle) as well as a bet- 
ter knowledge of the stellar masses. Nevertheless, the 
secondary star remains unseen. Even though the pre- 
dicted angular separation of 7 Cep B (0'.'84 for 2007.0; 
0'.'99 for 2009.0) is not particularly challenging, the 8- 
magnitude brightness difference in the visual band rela- 
tive to the glaringly bright primary explains all negative 
results (e.g., the speckle interferometry attempts by Ma- 
son et al. 2001, as well as the imaging by Hatzes et al. 
2003). We expect the contrast to be much more favor- 



able in the near infrared (Am ^ 6.4 in K), and that 
this detection should not be very difficult at those wave- 
lengths with adaptive optics on a large telescope. Such 
measurements of the relative position would allow a dy- 
namical determination of the mass of both stars, free 
from assumptions. 
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TABLE 1 

Heliocentric radial velocity measurements for 7 Cep derived in this work 
AS well as others collected from the literature, all placed on the CfA 

reference frame. 
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TABLE 2 

Ground-based positional measurements of 7 Cep from transit circle and photographic 
programs, on the international celestial reference frame. 
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^Includes scale factors described in the text. 
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TABLE 3 

Global orbital solution for 7 Cep. 



Parameter Value 



Adjusted quantities from outer orbit (A+B) 

Pab (days) 24392 ± 522 

Pab (yr) 66.8 ± 1.4 

7 (kms-i) -42.958 ± 0.047 

Ka (kms-i) 1.925 ± 0.014 

EAB 0.4085 ± 0.0065 

UJA (deg) 160.96 ± 0.40 

Tab (HJD-2,400,000) 48479 ± 12 

Tab (yr) 1991.606 ± 0.032 

a'l (mas) 324.6 ± 8.4 

JAB (deg) 118.1 ± 1.2 

^AB (deg) 13.0 ± 2.4 

Adjusted quantities from inner orbit (Aa+Ab) 

Pa (days) 902.8 ± 3.5 

Pa (yr) 2.4717 ± 0.0096 

/CAa (ms-i) 27.1 ± 1.5 

EA 0.113 ± 0.058 

'-'Aa (deg) 63 ± 27 

Ta (HJD-2,400,000) 53146 ± 72 

Ta (yr) 2004.38 ± 0.20 

Other adjusted quantities 

ARVi (kms-i) [McDonald I]'' . . . -45.228 ± 0.035 

ARV2 (kms-i) [McDonald 11]=^ . . -45.424 ± 0.035 

ARVi (kms-i) [McDonald III]'' . -44.053 ± 0.035 

ARV4 (kms-i) [CFHT]'' -44.483 ± 0.035 

ARV5 (kms-i) [CfA]'' +1.13 ± 0.12 

Ao* (mas) +73.6 ± 7.5 

A(5 (mas) +160.1 + 3.9 

A^i* (mas yr-i) -16.0 ± 1.1 

A/is (mas yr-i) +21.91 ± 0.81 

Derived quantities 

R.A. (sec)'" 23''39™ 21.0050 ± 0.0023 

Dec. (arcsec)*" +77°37' 55.241 ± 0.004 

/X* (mas yr-i) -64.8 ± 1.1 

lis (mas yr-i) +149.09 ± 0.81 

TT (mas) 72.70 ± 0.39 

a^g (arcsec) 1.382 ± 0.047 

aAB (AU) 19.02 ± 0.64 

/(A/b) (Mq) 0.01371 ± 0.00049 

Mb (Mq)" 0.362 ± 0.022 

/(Mp) (10-9 M0) 1.83 ± 0.32 

Mp sin iA (Mjup)= 1.43 ± 0.13 

aAa-p (AU)'^''^ 1.94 ± 0.06 



^Offsets to be added to the corresponding data sets in order to 
place them on the GrifTln system. 
''Coordinates of the barycenter (ICRF, J2000, epoch 1991.25). 
^Assumes a primary mass of MAa ^ 1-18 ± 0.11 M0 (sec g}. 
^Relative scmimajor axis of the orbit of the substcllar companion. 
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TABLE 4 

Spectroscopic determinations of the effective temperature and metallicity for 7 Cep from 

THE literature. 
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Source 
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ISoinrad & Taylor f 19691 

[Bakos (1971) 

jaicbocki (1972) 

^^istaffson, Kiacrgaa^^^^^idereen (1S7^ 

Camp bell (1978) 

Lgmber^^^^^ ^^^) 

Gratton et al. (1982) .. 

l^ae^aar^^T ^^^^^82 *) 

[Gratton (1985TT 

[Brown ct al. ( 1989) 

'McWilliam (1990) 

[Luck & Challcncr (1995) 

^Mishcnina, Kutsenko & Musaev ^^9^ . . . . 

[Soubiran. Katz & C avrcl (1998T 

[Gray ct al. (2003) 

^Santos. Israclian^^fe^r ^224) 

[Pranchini ct al. (2004) . ~ 

'Puhrman n (2004) 

Affer et a l. (2005) 

Luck fc Hcitcr (20Q^) 

This paper 



4383"^ 



+0.05^^ 
+0.02 
+0.1: 



4421*^ 


—0.04 






(4828) 


—0.21 ± 


25 


3.3 


(4630) 


+0.04 ± 


15 


3.1 


4840 


+0.02 ± 


08 




5091 ± lOO'^ 


-0.05 ± 


18 


3.57 ± 0.46 


4825 ± 60 


-0.04 ± 


14 


2.77 ± 0.15 


(4790) 


+0.04 




3.1 




-0.06 ± 


12 


2.77 


(4720) 


-0.04 




3.1 


(4770) 


0.00 ± 


11 


3.27 ± 0.40 


(4650 ± 100) 


-0.02 ± 


10 


2.35 ± 0.25 


(4810 ± 100) 


-0.02 ± 


10 


3.00 ± 0.30 


4769 ± 86 


-0.01 ± 


16 


2.98 ± 0.28 


4761 ± 80 


+0.07 ± 


08 


3.21 


4916 ± 70 


+0.16 ± 


08 


3.36 ± 0.21 




-0.066 ± 


034 




4888 ± 80 


+0.18 ± 


08 


3.33 ± 0.10 


4935 ± 139 


+0.14 ± 


19 


3.63 ± 0.38 


5015 ± 100 


+0.26 ± 


11 


3.49 ± 0.10 


4800 ± 100 






3.1 ± 0.2 



Note. — When not rep orted in the o riginal publications, typical uncertainties for [Fe/H] have been assumed 
to be 0.1 dex (0.25 dex for lSoinrad & Ta ylor 1969), and uncertainties in the effective temperatures have been 
assumed to be 100 K. 

^Temperature estimates given in parentheses arc listed for completeness, but are photometric rather than 
spectroscopic, and arc not considered further. 

^Although these values are listed as effective temperatures in the catalog by 
ICavrel de Strobel. Soubiran & Ralit3 i200lD , they arc actually excitation temperatures. We do not use 
them here. 

^The original value reported is +0.27. However, examination of the iron abundances derived for 12 other 
stars in this study indicates the [Fe/H] values are systematically overestimated by approximately 0.22 dex. 
Correcting for this offset brings the estimate for 7 Cep more in line with the rest of the determinations. We 
adopt the revised value here. 

"^The hotter temperature derived in this study is a consequence of the use of old values of the oscillator 
strengths fsee lMishenina. Kutsenko & Musae VI1995I1 . We have elected not to use it here. 
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TABLE 5 

Photometric estimates of the effective temperature 
OF 7 Cep. 



Photometric system and index T^ff (K)'' 



Johnson (B - V) 4756 ± 53 

Stromgren {b - y) 4811 ± 76 

Vilnius (Y -V) 4753 ± 79 

Vilnius {V -S) 4741 ± 70 

Geneva {B2 - VI) 4772 ± 51 

Geneva (52 - G) 4746 ± 44 

Geneva {t = [B2 - G] - 0.39[B1 - B2]) . . . . 4729 ± 49 

Johnson-Cousins {V - Rc) 4696 ± 73 

Johnson-Cousins (V - Ic) 4783 ± 52 

Cousins {Rc - Ic) 4893 ± 93 

DDO C(42 - 45) 4672 ± 63 

DDO C(42 - 48) 4729 ± 54 

2MASS {V - J)"" 5032 ± 370 

2MASS {V - H)^ 4972 ± 196 

2MASS {V - K)^ 4886 ± 209 

Tycho {Bt -Vt) 4749 ± 83 

Tycho-2MASS (Vt - K)^^ 4876 ± 194 



^Based on the color/temperature calibrations by 



IRamirez* & Melende^ I2005D for giants, adopting [Fe/H] ^ 
+0.01 ± 0.05 and no reddening (sec J4]. 

'^Due to the brightness of 7 Cep the star was saturated in the 
2MASS measurements and yielded a large photometric error. 
This is reflected in the large temperature uncertainty. 
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Fig. 1. — Radial velocity measurements of 7 Cep A as a function of time, along with our fitted curve from the combined solu- 
tion. The center-of-mass velocity o f the system is indicated by the dotted line, and is on the reference frame of the velocities by 
IGriffin. Carguillat fc GinesteU 12003) . (a) Classical velocity measurements in the outer orbit. The wiggles in the curve correspond to 
the perturbation by the 2.47-yr substellar companion, (b) Close-up of the high-precision velocities, which are near periastron passage in 
the outer orbit (arrow). The error bars in this panel are smaller than the size of the points. 
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Fig. 2. — High-precision radial velocity measurements of 7 Cep shown as a function of phase in the inner orbit, along with our fitted curve 
from the combined solution. The motion in the outer orbit has b een subtracted. The center-of-mass v elocity of the system is indicated by 
the dotted line, and is on the reference frame of the velocities bv lGriffin. Carguillat fc Gineste^ 12003) . 
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Fig. 3. — Path of 7 Cep A on the sky resulting from the combined effects of proper motion, orbital motion (P = 66.8 yr), and parallactic 
motion (solid curve). The magnitude and direction of the annual proper motion is indicated by the arrow. The Hipparcos observations 
are shown as dots at their predicted locations, and do not represent the actual measurements, which are one-dimensional in nature (see 
text). The dotted curve shows the path the star would follow in the absence of orbital motion, starting at the epoch of the first Hipparcos 
measurement, indicated with the open circle. 
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Fig. 4. — Computed orbit of 7 Cep A around the center-of-mass of the binary (shown with a plus sign). The direction of motion 
(retrograde) is indicated by the arrow, and the dotted line represents the line of nodes. The Hipparcos observations are displayed with 
filled circles at their predicted locations, and are seen to bracket periastron passage (open circle labeled "P"). The perturbation due to the 
substellar companion of the primary is negligible on the scale of this figure. The relative orbit of the binary is simply a scaled-up version of 
the ellipse shown here, with a scale factor given by (Mas + M^)/M^ = 4.26 (yielding a semimajor axis a^g = 1{'382 ± 0"047; see Table 151. 
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Fig. 5. — Enlargement of Figure HI showing the individual Hipparcos observations. See text for an explanation of the graphical represen- 
tation of these one-dimensional measurements. 
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Fig. 6. Ground-based catalog positions of 7 Cep in Right Ascension and Declination, after subtracting the contribution from the proper 
motiou resulting from our fit between the date of each observation and the reference epoch 1991.25. The curve represents the combination 
of motion in the 66.8-yr binary orbit and the parallactic motion, as predicted from the solution. 
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Fig. 8. — C onstraint on the orbital period of the binary for two different data sets: the one in this paper (sohd Hne), and the one used 
bv lGrittin. C arguillat & Gincstct 12002) (dashed hne). In each case the period has been fixed over a fine grid of values and the remaining 
elements were solved for in the usual manne r . The run of the reduced values for the present solution shows that of the two periods 
allowed by the lGriffin. Carguillat &: Gineste3 120021) fit, the 66-yr value is the correct one. 
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Fig. 9. — Evolutionary tracks from the calculations bv IYi et al.l 1200 1') and'Dema raue et alJ 120041) in the absolute visual magnitude vs. 
effective temperature plane. The metallicity adopted is [Fe/H] = +0.01, the weighted average of all spectroscopic determinations. Masses 
are labeled in solar units, and the dot with the shaded error box represents the measurements for 7 Cep. An enlargement is shown in the 
inset. 
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Fig. 11. — Theoretical mass and age combinations that are consistent with the four measured properties of 7 Cep {T^g, [Fe/H], My, and 
R) within their errors. The best fit is for = I-ISIq'J^ Mq and an age of 6.6^Q y Gyr. The larger point sizes indicate a closer match. 
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C78 = Campbell (1978) 
G83 = Graton et al. (1982) 

S98 = Soubiran, Kalz & Cayrel (1998) 

G03 = Gray et al. (2003) 

S04 = Santos, Israelian & Mayor (2004) 

F04 = Fuhrmann (2004) 

AOS = AHer et al. (2005) 

LOS = Luck & Heiter (200S) 
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Fig. 12. — Effective temperature and metalliclty determinations for 7 Cep from only the studies that measured both quantities spectro- 
scopically (filled circles). There is an apparent correlation between [Fe/H] and T^ff. The values adopted in this paper are indicated with an 
asterisk. The arrows represent other spectroscopic abundance determinations that do not have a corresponding spectroscopic temperature 
measurement. 
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Fig. 13. — Confidence levels of orbital fits in the iA-^A space of parameters, describing the wobble of 7 Cep A on the plane of the sky in 
response to the pull from the substellar companion. We use these fits to place an upper limit on the mass Mp of the substellar companion. 
The light gray area represents solutions for fixed values of these two parameters that can be ruled out only at the ~68% confidence level 
(Iit) or less. The plus sign corresponds to the best fit, which however does not give a statistically significant result (see text). The middle 
shade of gray corresponds to fits ruled out at confidence levels between Icr and 2a, and yields an upper limit on Mp of 13.3 Mjup. The dark 
gray area corresponds to fits ruled out at confidence levels between 2a and 3a. The outer edge of this region yields a minimum inclination 
angle of 4?9, and an upper limit on Mp of 16.9 Mjup. 
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APPENDIX 



ZERO-POINT CORRECTIONS TO THE RADIAL VELOCITIES OF 7 Cep FROM THE 

LITERATURE 



The historical sources containing radial velocity measurements of 7 Cep typically include other stars observed either 
as standards or for other purposes. The likelihood that many of those stars have been observed multiple times at 
the CfA is fairly high given that the spectroscopic database at CfA contains tens of thousands of stars and about a 
quarter of a million spectra to date. This common ground enables us to place the measurements of each of the sources 
on the CfA velocity system. In each case we selected all stars with no obvious signs of velocity variation that have 
been observed at least 3 times at the CfA. Radial velocities were derived from the available CfA spectra in the same 
way as those for 7 Cep, by cross-correlation usi ng synthetic templates based on model atmospheres by R. L. Kurucz 
f see iNordstrom et al.lll994 iLatham et al.ll2002D . The optimal template for each star was determined from grids of 
cross-correlations against a large number of synthetic spectra over broad ran ges in th e template parameters (mainly 
the effective temperature and rotational velocity), in the manner described bv lTorres. Neuhauser & Guenthcr (2002). 
Solar metallicity was assumed throughout. Many of the stars are giants but there are some dwarfs as well, so the 
optimal surface gravity for the template in each case was determined by repeating the procedure above for a range of 
values of \ogg, and selecting the one giving the highest correlation averaged over all exposures of the star. The radial 
velocities derived with these templates were then compared with those from each literature source. 

Some of these sources have relatively few stars that can be used as standards, and rejecting objects that have not 
been observed at CfA leads to the loss of potentially useful comparison stars in some cases, which can compromise the 
determination of the offset. Those stars can still be used so long as they are included in another of the data sets, which 
then provides the link to the CfA system. Thus, instead of separately comparing each source with CfA to determine 
the corresponding velocity offset, as might commonly be done, we have followed a procedure by which we determine 
the velocity offsets of all sources simultaneously by minimizing the scatter of the velocities for all standard stars taken 
together. In this way any star that is included in at least two of the data sets (whether or not one of them is CfA) 
can be used to strengthen the solution. The quantity we seek to minimize is 



where the sums are performed over all data sets (z = 1, . . . , -/Vgets)) stars in each data set (j = 1, . . . , iV^ star), and all 
observations of each star {k — 1, . . . , Ay^obs)- The quantity aijk represents the uncertainty of each observation. The 
mean radial velocity for each star, RVj, is a function the adjustable parameters (offsets A^) given by 



and changes as the iterations proceed. Since the offsets are computed relative to CfA (defined here as the first data 
set), Ai = 0. 

Table IX6l presents the results for each data set from our least-squares solution. We list the derived offset along with 
its uncertainty, the number of standard stars in each group, the number of observations of 7 Cep, and the interval 
of those observations. With a few exceptions the total number of standard star observations used in each data set 
is typically a few dozen, while the overall number of CfA observations used for those same standards is ~3300. The 
offsets were added with their corresponding sign to the individual velocities of 7 Cep in each data set to place them 
on the CfA system. 
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TABLE A6 

Radial velocity offsets applied to the literature sources containing 

MEASUREMENTS OF 7 CeP, TO BRING THEM ONTO THE CfA SYSTEM. 



Source 



Frost fc Adams (1903^ 

^cloDolskv (1904^ 

^liDhcr (1905) 

Kiistncr { 1908) 

Tibt (1973) 

Harper ( 1934) 

Joulon (1957) 

^nowden & Young (2005) . . . 

Kiacrgaard et al. (1981) 

TifA 



Oflfset A 


Standard 


RVs for 


Time span 


(kms-i) 


stars 


7 Ccp 


(yr) 


-1.33 ± 0.47 


12 


3 


1902.8-1903.3 


-0.48 ± 0.93 


5 


4 


1903.7 


-0.97 ± 0.58 


9 


3 


1905.8 


-1.38 ± 0.70 


12 


5 


1905.9-1907.8 


+0.46 ± 0.23 


14 


3 


1916.0-1917.8 


+1.82 ± 0.22 


21 


1 


1921.9 


-1.50 ± 0.28 


6 


1 


1955.0 


+0.75 ± 0.16 


14 


13 


1972.5-1974.7 


+0.11 ± 0.48 


14 


1 


1977.7 


0.00 


71 


3 


2001.5-2004.9 



